### Frequentist logit model estimation w/o interaction ##

. logit no_shop cdc_m pres_m state_m expert_m health_frame shelter jobloss white education gop gender ideology_rs

### Frequentist logit model w/ interaction terms ###

. logit no_shop i.cdc_m i.pres_m i.state_m i.expert_m i.health_frame cdc_m##health_frame pres_m##health_frame state_m##health_frame expert_m##health_frame i.gender i.education i.white i.gop i.shelter i.jobloss ideology_rs

. margins

. margins health_frame#expert_m health_frame##pres_m

. margins health_frame#pres_m

. margins health_frame#expert_m

. margins, dydx(health_frame expert_m pres_m)

### Bayesian model estimation w/o interaction terms ###

. bayes, block({no_shop: cdc_m}) block({no_shop: pres_m}) block({no_shop: state_m}) block({no_shop: expert_m}) block({no_shop: health_frame}) block({no_shop: shelter}) block({no_shop: jobloss}) block({no_shop: white}) block({no_shop: education}) block({no_shop: gop}) block({no_shop: ideology_rs}) block({no_shop: gender})  mcmcsize(200000) burnin(20000) saving(simdata1) rseed(32306) prior({no_shop:}, uniform(-10,10)): logit no_shop cdc_m pres_m state_m expert_m health_frame shelter jobloss white education gop gender ideology_rs 
  
### Bayesian model estimation w/ interaction terms ###

. bayes, block({no_shop: cdc_m}) block({no_shop: pres_m}) block({no_shop: state_m}) block({no_shop: expert_m}) block({no_shop: health_frame}) block({no_shop: shelter}) block({no_shop: jobloss}) block({no_shop: cdc_frame_h}) block({no_shop: pres_frame_h}) block({no_shop: state_frame_h}) block({no_shop: expert_frame_h}) block({no_shop: white}) block({no_shop: education}) block({no_shop: gop}) block({no_shop: ideology_rs}) block({no_shop: gender})  mcmcsize(200000) burnin(20000) saving(simdata2) rseed(32306) prior({no_shop:}, uniform(-10,10)): logit no_shop cdc_m pres_m state_m expert_m health_frame shelter jobloss cdc_frame_h pres_frame_h state_frame_h expert_frame_h white education gender gop ideology_rs

### Bayesian Interval hypothesis tests ###

. bayestest interval {no_shop: health_frame}+{no_shop:cdc_frame_h}, lower(0)

. bayestest interval {no_shop: health_frame}+{no_shop:pres_frame_h}, lower(0)

. bayestest interval {no_shop: health_frame}+{no_shop:state_frame_h}, lower(0)

. bayestest interval {no_shop: health_frame}+{no_shop:expert_frame_h}, lower(0)